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Two-dimensional loop-erased random walks (LERWs) are random planar curves whose 
scaling limit is known to be a Schramm-Loewner evolution SLE K with parameter k = 2. In 
this note, some properties of an SLE K trace on doubly-connected domains are studied and 
a connection to passive scalar diffusion in a Burgers flow is emphasised. In particular, the 
endpoint probability distribution and winding probabilities for SLE2 on a cylinder, starting 
from one boundary component and stopped when hitting the other, are found. A relation of 
the result to conditioned one-dimensional Brownian motion is pointed out. Moreover, this 
result permits to study the statistics of the winding number for SLE2 with fixed cndpoints. 
A solution for the endpoint distribution of SLE4 on the cylinder is obtained and a relation 
to reflected Brownian motion pointed out. 
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I. INTRODUCTION 

The Schramm-Loewner Evolution SLE K is a one-parameter family of conformally invariant mea- 
sures on non-intersecting planar curves. In many cases, interfaces in two-dimensional models of 
statistical mechanics at criticality are conjectured to be described by SLE. For example, the scaling 
limit of the planar loop-erased random walk (LERW) is known to be SLE2, interfaces in the 2D 
Ising model SLE3 and critical percolation hulls to be described by SLEg. While these systems 
can be also studied using traditional methods of theoretical physics like (boundary) conformal 
field theory and Coulomb gas methods, Schramm has given a novel and rigourous approach via 
conformally invariant stochastic growth processes. 

LERWs were the starting point and motivation for the development of SLE in Schramm's seminal 
paper [l| . Today their scaling limits and relations to SLE are well understood in simply-connected 
domains, where conformal invariance was shown by Lawler, Schramm and Werner [2j. Properties 
in multiply-connected domains are more involved but much recent progress was obtained, e.g. in 
a series of works by Dapeng Zhan a a a. In this paper, we mainly study properties of planar 
LERWS/SLE2 in doubly-connected domains, emerging from one boundary component and condi- 
tioned to finish when hitting the second boundary component. In contrast to the simply-connected 
case, curves in doubly-connected domains may wind non-trivially around various boundary compo- 
nents. Here we concentrate on these winding properties of loop-erased random walks on cylinders 
in two cases: (1) with one endpoint fixed on one boundary component, and the other endpoint free 
on the other boundary component, (2) with both endpoints fixed on different boundary compo- 
nents. We establish these results by solving a differential equation arising from SLE K with k = 2 
on the cylinder. Another motivation to study the cylinder is that numerical simulations of critical 
systems are often performed in a cylinder geometry, see e.g. 
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The paper is organised as follows. In section 2 we give a brief description of SLE. In particular, 
we outline two variants of SLE in simply-connected domains which prove to be limits (in some 
sense) of the doubly-connected case. After this we discuss SLE on cylinders, being prototypes 
of doubly-connected domains. Section 3 provides a detailed discussion on endpoint probabilities 
and winding properties of SLE K on cylinders. Inspired from analogies with diffusion-advection in 
(1 + l)-dimensional Burgers flows, we present analytical results for the case k = 2, study various 
limits and point out a relation between the winding of LERWs and conditioned one-dimensional 
Brownian motion. Finally, conditioning walks to exit via a given point, we derive the winding 
statistics for LERW with fixed endpoints on the cylinder. In section 4, somewhat aside the main 
topic, we provide the endpoint distribution in the case k = 4 and point out a relation to reflected 
two-dimensional Brownian motion. We conclude in section 5. The appendix contains computations 
of exit probabilities for Brownian motion and random walks on the lattice, in relation to distribution 
of exit points for SLE2 and their finite-size corrections. It also contains a partial (periodised) result 
about bulk left-passage probabilities with respect to a point in the covering space of the cylinder, 
as well as a path-integral version of the endpoint distribution problem. 



II. BASIC NOTIONS OF SLE 

In this section we give an informal description of SLE. For detailed accounts we refer the reader to 
the many existing excellent review articles on this subject (see e.g. 0,0, @])- 

Stochastic Loewner evolutions (SLE) describe the growth of non-intersecting random planar curves 
or hulls emerging from the boundary of planar domains B by means of stochastic differential 
equations. A crucial feature of SLE is conformal invariance: if two domains B and B' are related 
by a conformal mapping g(z), then the image of SLE in B is SLE in B'. Another key property is 
the so-called domain Markov property of SLE. Conditioning on the curve/hull jt grown up to time 
t in B, the remainder of the curve is just SLE in the slit domain H>t = D\jf In fact, these two 
properties are (almost) imply that the growing hulls a characterised by a single positive parameter 

K. 

For simply-connected domains, conformal invariance and the Riemann mapping theorem naturally 
lead to a description of the curves in terms of conformal mappings gt{z) which uniformise Bt to B, 
and map the tip of 74 to some point on the boundary of B. gt(z) is solution of a differential equation 
in t with initial condition gt=o(z) = z, called the Loewner equation. SLE in simply-connected 
domains comes in different flavours, according to which the Loewner differential equation varies. 
In section Hi A\ we shall concentrate on so-called radial and dipolar SLEs. 

Compared to simply-connected domains, SLE in multiply-connected domains, proves to be more 
involved because no equivalent of the Riemann mapping theorem is at hand [3|, 14, LLCJ, 111]- In this 
article, we concentrate on doubly-connected domains B with the topology of an annulus. Their 
boundary is given by two disjoint simple curves <9iB and c^B. A classical theorem of complex 
analysis states that each such a domain may be mapped onto the annulus A p = {z £ C | e _p < 
\z\ < 1} with some p > 0, such that for example diO is mapped to \z\ = 1 and c^B to \z\ = e~ p . 
Conformal invariance of SLE legitimates the choice of A p as a reference domain. However, we shall 
rather work with a cylinder T p = {z = x + iy G C | — tt <x<ir,0<y<p}, identifying x = —tt 
and x = tt. The cylinder T p may always be mapped back onto A p via the conformal mapping 
w = e lz . The height of the cylinder p > is called the modulus, and the boundary components are 
given by d\T p = [—it, tt] and c^Tp = [—tt, tt] + ip which we shall refer to as "lower" and "upper 
boundary" respectively. We shall study hulls jt growing from z = G d\T p to 82 T p as shown 
on figure QJi and b. In section III Bt we analyse the Loewner equation for gt(z) mapping T p \yt 
to some cylinder T q with (in general) different modulus q 7^ p. A mapping back to T p does not 



exist since one may show that there is no conformal mapping sending a cylinder of modulus p to 
a cylinder of modulus q 7^ p 12|]. Moreover, we show how the doubly-connected version of SLE on 
T p interpolates between radial and dipolar SLE. 



A. SLE in simply-connected domains 

Before embarking into considerations on general values for the modulus p we shall consider limiting 
cases of semi-infinite cylinders p — > +00 and very thin cylinders p —* + in order to establish a 
connection with well-known variants of SLE in simply-connected domains. 

p — > +00 and radial SLE. Radial SLE starts from a boundary point xo G <9B to a bulk point 
Xoo E D of some simply-connected planar domain B. Its Loewner equation is most conveniently 
written for the unit disc B = U from x$ = 1 to Xqo = where 

T7 = -9t[.z)—— 1) 

at gt{z)-e l & 

where £j is some real- valued function of t > and gt(0) = (the point Xqo remains unchanged 
during time evolution). The time parametrisation has been chosen so that g' t (0) = e . The real- 
valued function £ t is related to the image of the tip Tt of the growing trace jt by gtijt) = e ^'- I n 
fact, conformal invariance and the domain Markov property imply that £j = \J~kBt + at is simple 
Brownian motion with a diffusion constant k and a linear drift at. If we ask for rotation invariance 
as it occurs in physical systems (and throughout this paper), the drift vanishes a = 0. For later 
considerations, we outline another useful version of radial SLE. In fact, the unit disc U may be 
sent to the semi-infinite cylinder T ro via the conformal mapping w = — ilogz (here the principal 
branch of the logarithm has been chosen). In this geometry, the Loewner equation for the new 
uniformising map gt( z ) is given by 

-^r = cot { — 2 — J' &=^ B f ( 2 ) 

obtained from ([T]) using gt(z) = exp[ig t (w(z))]. 

Since we aim at studying winding probabilities, let us recall how to estimate the distribution of 
the winding angle 6t for large t [H, 0]. For the unit disc U as t > 1, the tip of the radial SLE 
trace Tt approaches z = 0. We may therefore write gt(Tt) = e 1 ^ ~ e*rj, by using the Taylor series 
expansion gt(z) = e t z + . . . near z = 0, from what follows Tt ~ e _<+ ^*. Hence the winding angle 
9t = argTf = £t turns out to be a Gaussian random variable with zero mean and variance Kt. For 
the infinite cylinder the winding angle corresponds to the real part of the tip allowed to vary 
continuously on R, and its probability distribution is approximately given by 

1 fx 2 



p — > + and dipolar SLE. Another version of SLE is obtained by considering planar curves in 
a simply-connected domain B from a boundary point xq G <9B to a side arc S C <9B [l3| not 
containing xq. Let us choose the geometry of an infinite strip § p = {z G C | < Im z < p} of height 
p, xq = and S 1 = R + ip. Then Loewner's equation is given by 



%^ = ^coth 
at p 



n(gt(z) - & 

2j9 



(4) 



gt(z) maps the tip of the trace 74 to £t. Again, because of conformal invariance and the domain 
Markov property we have = \J~kBi. The trace 74 hits the boundary arc M + ip at some (random) 
point 2; + ip as t — > +00. This version of SLE is related to the cylinder problem on T p in the limit 
p — > + . Loosely spoken, if p is very small the SLE trace will not feel the periodicity of the cylinder 
but only hit the upper boundary close to the point ip on the upper boundary (this statement can 
be made more precise, see below). The equivalent to the winding angle distribution for radial SLE 
is the probability distribution of this exit point for dipolar SLE, which we shall denote \d(p, x). 
For any k it is given by [13| 

_ v^r(l/2 + 2/K) 1 
MM 2pT(2/ K ) [cosh(7rx/2p)] 4 / K [) 

B. Stochastic Loewner evolution on doubly-connected domains 

Let us now turn to the case of general modulus p. The Loewner equation for the uniformising map 
g t (z) of SLE K on T p is given by 0] 

^M = H(p-t,g t (z)-V^B t ) (6) 
with standard Brownian motion Bt and the function 

n 00 „ . 

tti \ t ( z -i \ z sr^ 2smz 

H(p,z)= hm > cot ikp = cot — h > : — ; 7 

' n^oo ^ V2 V 2 ^ cosh 2/cp - cos z w 

k=—n k=l 

The process is defined between i = (50 (z) = z) where the trace starts from point z = on the 
lower boundary up to t = p where the SLE trace hits the upper boundary at some random point 
7 P with Im7p = p. The time parametrisation of the curves is chosen in such a way that gt(z) maps 
T p \7j to T q with q = p — t. One possible way to derive © and ([7]) is by considering an infinitesimal 
hull, and using invariance properties of T p , namely translation invariance z z + a, aSR and 
reflection invariance z 1— ► ip — z (see where this has been explicitly carried out for an annulus 
geometry) . 

Let us describe the properties of the function defined in ([7|). H(p, z) is an odd function in z with 
period 2ir. Furthermore it is quasi-periodic in the imaginary direction: H(p, z + 2ip) = H(p, z) — 2i. 
Moreover, notice the interesting property H{p,z) = (iirH(ir 2 /p,mz/p) — z)/p. It may be used in 
order to derive an alternative series expansion 

For large positive p, we have H(p — t,z) = cot(z/2) + . . . as it can be seen from ([7|). Therefore, 
as we let p — > +00, the Loewner equation Q reduces to the one for radial SLE, defined in ([2]). 
The dipolar limit p — > + is more involved. Consider the function ht(z) = gt{z) — yfUBt mapping 
Tp\jt to T p _ t and the tip r t of the curve to 0. Rescaling k s (z) = ir ht{z) /{p — t) maps the problem 
to a cylinder of fixed height ir, and width 2ir 2 /(p — t). The new mapping k s (z) depends on a 
time s = tt 2 /(p — t), < t < p. Performing the time change leads to a stochastic differential 
equation dk s (z) = iH(s,ink s (z))ds — y/KdB s where B s is another standard Brownian motion, 
obtained from the time-change formula. Now, as p — > + we have s — > +00 so that we may use 
H(s, z) ~ cot z/2 + . . . . Finally for K s (z) = k s {z) + t/kBs, we find the dipolar Loewner equation 
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Figure 1: (a) Illustration for SLE on an annulus A p from 1 to the target circle \z\ — e~ p . (b) SLE on the 
cylinder and the covering space: the shaded region represents the cylinder, with sides x = —n and x = n 
identified. We obtain the covering space by periodic extension along the real axis. The endpoints on the 
cylinder and the covering space differ by 2?tN where N denotes the winding number of the trace. The 
interfaces are samples of loop-erased random walks. 

in the limit s — > +oo. 

The covering space. In the sequel it shall be convenient to lift the processes to the covering space 
of the cylinder T p . In fact, this is simply obtained by allowing Re7t to vary continuously on R. 
An illustration is given on figure [TK It is easy to see that the difference between the endpoints of 
a given trace on the cylinder and its analogue on the covering space is always given by 2irN with 
an integer N. We shall call N the winding number of the SLE trace. For the covering space, the 
Loewner equation ([6]) still holds. However, probabilities and probability densities computed on the 
covering space have different boundary conditions when compared to the periodic case of cylinders 

a. 

III. ENDPOINT PROBABILITIES FOR SLE ON THE CYLINDER 

In this section, we characterise the endpoint distributions for SLE on cylinders of height p which we 
denote A(p, x), —ir < x < ir. To be precise, A(p, x)dx corresponds to the probability that the trace 
hits the upper boundary within Re 7 P G [x, x + dx] for some infinitesimal dx. In order to study the 
winding of the traces around the cylinder, we introduce an analogous object X(p, x), x £ R for the 
covering space. Results on the cylinder are obtained upon identification of points differing by lim: 

oo 

A(p,x)= V] X(p, x + 2im) (10) 

n=— oo 

A. Endpoint distributions for general n 

In order to compute A(p, x) consider the SLE trace 7dt evolving for an infinitesimal time dt on 
the cylinder T p . The uniformising map gdt(z) maps the remainder of the curve 7 p \7dt onto SLE 
in a cylinder T p _^ t starting at y^di^o- Under this mapping gdt( z ) the point x + ip on the upper 



boundary of T p is sent to some point x' + i(p — dt) on the upper boundary of T p _d t , with x' = 
x + (H(p, x + ip) + i)dt, see ([6]). Using the Markov property and conformal invariance of the SLE 
measure we obtain, upon averaging over the infinitesimal Brownian increment dBo, the equation 

A(p, x)dx = Ed_B [A(p — dt, x — yJHdBo)dx'] (11) 

An expansion to first order in dt, using Ed_B [di?o] = and Eds [di?o] = dt, leads to the Fokker- 
Planck equation 

dMp.x) d . , ... Kd 2 A(p,x) . . 

where we have introduced the drift function (compare with 0) 

AT 

7T 7t(j; + 2-7rn) x 
v(p,x) = H(p,x + ip) + i = — hm > tanh (13) 

P N -*°°r^N 2P P 

Like H(p, z) the drift function v(p, z) is an odd and periodic function with period 2ir. Furthermore, 
it is quasiperiodic in imaginary direction: v(p,z + 2ip) = v(p,z) — 2i. Notice that X(p, x) is a 
solution of (]12p as well, however with vanishing boundary conditions for x — > ±oo. Leaving aside 
the identification of the points x and x + 2tt, we ask for the probability 

/X 
dy\(p,y) (14) 

that the SLE traces hits the upper boundary to the left of a given point x + ip on the covering 
space. It is solution of the convection-diffusion equation 

du{p,x) _ , duj^x) K d 2 uj(p,x) 

~dp~ " V W> x )—&— + 2 dx 2 ' (15) 

with boundary conditions lim p _ >0 + u(p, x) = @(x), lim x .^ +00 u(p, x) = 1 and lim^-^-oo u(p, x) = 0. 
Here Q(x) denotes the Heaviside function which is 1 for x > 0, and otherwise. It is a remarkable 
property of the function v(p, x) to be a solution of the (1 + l)-dimensional Burgers equation 

dv(p,x) _ dv(p,x) | d 2 v(p,x) ^ 

dp ' dx dx 2 ' 



where p plays the role of time and viscosity is set to unity [4|. In fact, this equation even holds 
for complex x. Furthermore is can be shown that H(p, z) obeys the same equation by using the 
relationship between H(p,z) and v(p,z). Hence we interpret (]15p as the evolution equation of a 
passive scalar convected and diffusing in the Burgers flow — v(p,x), also described by a Langevin 
equation: 

dx p = —v(p, x p )dp + \fhidBp (17) 

The problem of finding SLE probabilities for arbitrary k can thus be reformulated as finding the 
density of a passive scalar convected by a ( par ticular) one-dimensional Burgers flow (along x). In its 



general formulation it is a classic problem [14], [l5|, [l6|, LL7|, [Ha, LL9(] characterised by the dimensionless 



Prandl number P [20], with P = 2/k for our problem. 

By means of the well-known Cole-Hopf transformation 2l|, [22| we rewrite 



v(p,x) = 2— \ne(p,x), 



(18) 



where e(p, x) satisfies the diffusion equation de(p,x)/dp = d 2 e{p, x)/dx 2 with solution 



e(p,x)= (-l) n exp (mx-n 2 p) = ^ exp — ( 19 ) 

n=— oo " n=— oo ^ P / 

corresponding to an initial condition e(0, x) = 2ir Y^=-oo S(x — Tr—2-Kn). In the language of Burgers 
equation this corresponds to an initial condition at p = + with a periodic set of shocks at positions 
xft = 2im. As the "time" p increases the shocks are broadened by the viscosity and acquire a width 
proportional to p, leading to the smooth periodic function (|13p . Finding the solution of (| 12|) for 
general k is a complicated time-dependent problem, and only results in special cases are obtained 
below (see appendix [D] for a path-integral formulation). 

In the limit of large modulus p S> 1 the drift term in (|12p vanishes like v(p, x) ~ 4e~ p sinx. 
Therefore we recover a simple diffusion equation for the endpoint distribution X(p, x) which, in 
terms of the scaling variable x / y/p, converges to the Gaussian result for radial SLE Q . However if 
one does not rescale x, there are non-trivial corrections and one expects instead that the cumulants 
converge at large p to k dependent constants, (x 2n ) c — np5 n ,2 — > -^2n,K> as shown explicitly below 
for k = 2. Going from T p to an annulus with radii R\ > R 2 = R\e~ p via the conformal mapping 
w = R\e lz (see figure [Th^ for an illustration) we obtain the winding angle distribution 6 = x for 
SLE starting from the outer boundary and stopped when first hitting the inner boundary (or the 
reverse if we admit reversibility of annulus SLE) . We expect its cumulants of the winding angle to 
have the form 

(0 2n ) c = Kln(R 1 /R 2 )5 n , 2 + X 2n , K (20) 

in the limit of large R\/R 2 . This differs from the standard result, also obtained from Coulomb 
gas methods 0, E3] or via the argument given for radial SLE, namely (9 2 ) = (4/g) ln(i?i/i?2) 
with 4/g = k (see (J3j) ) , by universal constants. One can compare to the corresponding cumulative 
distribution of winding angle for the standard two-dimensional Brownian motion with absorbing 
conditions on the boundary of the annulus (it can be obtained from dipolar SLE with k = 2, the 
direct derivation being recalled in IA lj) 

Ub ^ = 1 + e -j/ln( Rl /R 2 ) ( 21 ) 

which yields much larger cumulants (9 2n ) c ~ (l n (^i/^2)) 2n with prefactors given below in the 
context of dipolar SLE. 

As p — > + a simple argument allows to recover the dipolar SLE hitting probability ([5]) from the 
form of Burgers shocks, which are well separated in that limit. Let us suppose that we may restrict 
(fT3l) to a single shock, v s (p, x) = —x/p + ir/pt&Yih.(irx/2p), i.e. take an initial condition e(0, x) = 
2tt(5{x + t:)+5(x — 'k)). This gives a potential u(p, x) = 21oge(p,x) = — x 2 /2p+2 log cosh(7rx/2p) + 
const, which in the scaling region x ~ p strongly confines the diffusing particle described by (|12|) 
near the origin. One easily checks that the (adiabatic) Gibbs measure A a diab(^ 5 p) ~ e -( 2 / K ) u (p> x ) = 
C(p)e{x,p)~ A / K = Xd(x,p) satisfies (fT2"j) in the region x ~ p (the term d\(p, x)/dp on the l.h.s. being 
subdominant). One can be more precise and rescale £ = x/p. The new probability u>(p, £) = uj(p, x) 
is solution of the differential equation 

where v(p,£) = p(v(p,x) + £). Inserting the single shock profile one finds v s (p, £) = it tanh(7r£/2), 
i.e. a stationary velocity field. Hence we may convince ourselves that the stationary solution of 
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Figure 2: Illustration of loop-erasure on the cylinder for a square lattice, (a) A Brownian excursion from 
z = to the upper boundary y — p. (b) Loop-erasure of the same walk. Loops like (14, 15, . . . , 27), winding 
around the cylinder for several times, may be erased. 



(|22p is d x uj s tat = Arf, i.e. the convergence to the dipolar result © includes also (non exponential) 
subdominant terms in the scaling region x ~ p. Interesting exact solutions of (|22l ) for all p are 
easily obtained from a Laplace transform in r = 1/p or using methods as in [20], but do not 
obviously relate to SLE. 



B. Results for k = 2 and the winding of loop-erased random walks 



The Cole-Hopf transformation (I18p will lead to interesting simplifications allowing to treat explicitly 
k = 2. This case corresponds to the scaling limit of the LERW on the cylinder. First we recall its 
definition on the lattice |9( . Consider a rectangular lattice domain of lattice mesh 5 embedded into 
the cylinder domain T p . Start a simple random walk W from to Rez = y = p on with absorbing 
boundary conditions for y = (see figure^). Chronological loop erasure from W yields a simple 
path 7 on the lattice domain which is the LERW (see figured). Notice that the endpoints of W 
and 7 always coincide. However, the random walk may wind several times around the cylinder 
whereas its loop erasure does not. This makes it a non-trivial exercise to compute the winding 
properties of the latter (i.e. its endpoint distribution on the covering space). As the lattice mesh 5 
goes to 0, the paths W converge to two-dimensional Brownian excursions on T p whereas the paths 
7 converge to SLE2 on T p 

Since the endpoints of W and 7 coincide on the cylinder T p , we may find the endpoint distribution 
for SLE2 by a simple computation on Brownian excursion on the cylinder, starting from z = and 
stopped as soon as it hits altitude y = p (see lAip . The result is 

A 2(P.*) = h E ~, , , , „ PS) 



(„»„ (2^=)))- 



An equivalent result for random walks on a finite cylindrical lattice domain is presented in appendix 
IB! Indeed, A2(p,x) solves (I12D for k = 2 with periodic boundary conditions and may be related 
to the drift function v(p,x) via Aa(p,x) = (pdv(p,x)/dx + 1)/2tt [4|]. Moreover, it corresponds to 
the periodised endpoint distribution (0) for dipolar SLE with k = 2. 



In order to compute the winding probabilities and to find the distribution on the covering space 
X(p,x), we solve (I15D with boundary conditions as indicated above. The idea is to write u(p,x) = 
a(p,x)^(p,x) and impose that terms proportional to dip(p,x)/dx cancel, which, for arbitrary k, 
yields a(p, x) oc e(p, x)~ 2 l K ^2$\. It turns out that the function ip(p,x) is solution to 

dijj(p,x) _ (2 - k\ dv(p,x)^^ ^ + K d 2 j){p,x) ^ 



dp \ 2k J dx ' 2 dx 2 

Hence for k = 2 we have a simple diffusion equation dip(p, x)/dp = d 2 ip(p, x)/dx 2 with initial 
condition 

oo 

ip(0,x) = e(0,x)u)(0,x) = 2vr ^ 8{x - tt - 2irn)Q{x). (25) 

n=— oo 

In this case the solution therefore is immediate: 

/OO J I — 

_^ =e ~(-y?/^ i){ ^ y)= Py e ~(— (2n+l))V4 P . (26) 
-oo V4vrp V P ^ 

Thus, for k = 2 the endpoint probability reads 

, s E^oexp(-( 3 :-^(2n + l)) 2 /4p) E^=o exp (^n(x - ^(n + l))/p) 
2 E^-ooexp(-(x-vr(2n + l))V4p) E^=-oo exp (vrn(x - vr(n + l))/p) ' 1 j 

Note that w(p, x) = ip(p,x)/(ijj(p,x) + ^>(p, — x)). The endpoint distribution is given by 

(x—TT—2-Kn) 2 + (x+7r+27rm) 2 ' 



E^ m =o( m + n + 1 ) ex P , , 

A 2 (P, X) = ^ " ( 28 ) 

{pe{p,x)/Try 

This result is exact and allows to analyse for n = 2 how the SLE trace wraps around the cylinder. 
An illustration of the distributions is given in figure El Some remarks are at order. First of 
all, we have checked through tedious calculation that (f28j) inserted in (fTUj) correctly reproduces 
(I23|) . Second, the very particular form of (|27p . given by a ratio of propagators for one-dimensional 
Brownian motion, suggests that it might be obtained as a conditional probability for a simple 
one-dimensional process. Indeed, such a relation can be established (see section UlI Cj) . Moreover, 
notice an interesting relation to (boundary) conformal field theory. Setting q = exp(—2i: 2 /p), we 
see that the different terms in (|27|) contain q h hn+i where h\^ n = dim^i^ = nin — l)/2 are the 
scaling dimensions of the operators ipi n , creating n—1 curves at the boundary, for the 0(N = —2) 



model 24|. Finally, one may hope that the simplifications in the case k = 2 for the endpoint 
probabilities may be extended to bulk probabilities of the cylinder and the covering space, a case 
which we treat (partially) in appendix [Q 

The closed forms (127)) and (j28H converge well in the dipolar limit p — » + , however are not suitable 
for computations in the radial limit p — * +oo. Therefore, we study the moment generating function 
{e~^ x ) here below, and shall give an exact expression adapted to both limits. To this end, it will 
be useful to consider L02(p,z) = ip(p, z)/e(p, z) for complex z. From ([27]) we see that it is periodic 



with period 2ip. Moreover, since the denominator can be written as an infinite product 25J] 



e(p, z) oc f[ (l - e ^-(2fc-i) P \ ^ _ e - iz -(2k-i)p\ (29) 



k=l 



X 2 (p,x), A 2 (p,x) 
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Figure 3: The endpoint distributions A.2(p, x) and \2{p, x) for SLE2 in the case p = n. 



(up to a p-dependent prefactor), we conclude that the function 0J2 (p,z) has simple poles in the 
complex plane for z^/ = (2k — l)ip + 2ir£, k, £ £ Z. In order to compute the generating function, 
we write 



fi > 0. 



(30) 



Deformation of the integration contour in the complex plane from R to 1 + 2ip leads to two 
distinct contributions. First, the integration along R + 2ip simply yields e~ 2tpil (e~ IMX } because of 
the periodicity of u>2(p,z). Second, we must take into account the simple poles at Zg = = 
ip + 2ir£, £ £ 7L and compute the corresponding residues. Solving for (e~^ x ), we find 



tix\ 



2ttiu \ -« 
7T. — > res e 

l&Z 



TTfJ, ^— v _2ttu^ 11 ■ , ( 9e(p, ip + 2tt£) 



sin pfi 



Y^e- 2nfie if>(p,ip + 2ir£) 



Oz 



(31) 



The sum can explicitly be computed because the derivative of e(p, z) at z = Z£ does not depend on 
t by periodicity. After a little algebra we obtain 



Pli 



1 



E^ (-l) £ smh((2l + l W )e- 



-w 2 e(e+i)/p 



sinp/i sinhvr^ YX=v(- l Y( 21 + l)e-^ +1 VP 
For large p it is convient to transform this expression into a dual series 

e^ 2 7r/i 1 YZv(~ l Y S H&£ + l)pn)e- pli - l+ ^ 



sinhvr^ sinpn EZoi^Y^ + l)e"P^ +1 ) 



(32) 



(33) 



by using identities for Jacobi theta functions [25| or Poisson's summation formula. Let us now 
study the radial and dipolar limit of our result. 



Limit p — > +00. In the radial limit, we may write (133]) as 



PA* 



7T/X e 

sinli7r/u 



;i+2e" 2p (l- cos 2pfi)) + 0(e" 



(34) 



In fact, it is possible to recover this result by approximation of e(p, x) with its series expansion (|19j) 
and Gaussian integration. The result turns out to be finite result for all [x since the probability 



decays as a Gaussian at large x, i.e. 1 — ^(p, x) = ^(p, — x) ~ exp(— (x + 7r) 2 /4p) as x — ► +00, 
up to a prefactor periodic in x. We take the logarithm and find a cumulant expansion 



00 ^ / „. \ 

log(e-^) = ^(-l) fc (x fc ) c |j = pm 2 + log ( -r?g— ) + 2e" 2 f (1 - cos 2p M ) + O (e"^) 



(35) 



The result shows that, as announced above, we do not recover a Gaussian distribution in the limit 
p — > +00, in contrast to the simple argument from radial SLE. It rather provides all the constants 
X2n,n=2 from the Taylor series expansion of the function log[7r/i/(sinh7iyz)] with respect to fj,. For 
instance, we find 

(x 2 ) c = 2p- vr 2 /3 + 8pV 2p + O (e~ 4p ) (36) 
(x 4 ) c = 2^/15 - 32pV 2p + O (e~ 4p ) (37) 

up to corrections of the form p n e~ 4p . 

Limit p —* + . In the dipolar limit, we expand the cumulant generating function, using (I32|) . as 



log( e -^) = log - 4sinh 2 (vr^)e-^ / p . (38) 

Hence the first two non-vanishing cumulants are given by 

(x 2 ) c = p 2 /3 - 8vr 2 e- 2 " 2 / p + O (e- 4w2 / p ) , (39) 

(x 4 ) c = 2//15 - 32vr 4 e- 2 " 2 / p + O ( e - 4 " 2 / p ) . (40) 

Since the first corrections to the dipolar result at small p originate from erasure of a single loop 
wrapping around the cylinder, it seems plausible that the probability for erasing such a loop 
behaves like exp(— 2ir 2 /p) in the limit p — > + . This factor has indeed the same magnitude as the 
probability that planar Brownian motion on T p winds by at least 2ir around the cylinder (we recall 
the winding distribution for Brownian motion in appendix I A ip . 

Notice that the dipolar limit presents a subtlety, which does not allow to conclude immediately 
from expansion of the probability u>2(p, x). For any fixed x such that — ir < x < tt we find estimate 
as p — ► + (a neighbourhood of size ~ p of the points ±7r must be excluded): 

u 2 (p, x) = TT A__ + tanh (™) e -»<>/P + o(e -- 2 / P) . (41) 

Up to exponential corrections, and uniformly in the interval ] — tt, tt[, we recover the dipolar result 
([5]) for k = 2. However, the approximation is not uniform beyond this interval what leads to 
complications when computing moments by using this approximation. 



C. Probabilistic argument for k = 2 

The simplifications for SLE2 seem surprising at first sight. Here we shall try to understand the 
underlying mechanism by giving a probabilistic argument for the simple solution of (I15p by (|27p 
in the case k = 2. 

Let us consider the stochastic process X t = g t (x + ip) — i{p — t) — \fnB t with initial condition 
Xq = x. It describes the motion of a given point x + ip on the upper boundary of the cylinder 
under the flow gt{z) and is solution of the stochastic differential equation 



dX t = v(p - t, X t )dt + v^dBf. 



(42) 



In particular, it has been shown that lim^- Xt = (2k + 1)tt with k £ TL (almost surely) |4|. We 
may relate the probability that the SLE trace hits the upper boundary of the covering space within 
(— oo, x] to Xt via 

tv(p,x) = P[Re7 p e (-00, x]) = P[Reg t (j p ) £ (-00, Reg t (x + ip)]] 

= P[V^B P e (-oo,Re 5p (x + ip)]} = P X [X P > 0], (43) 
t-+p 

where P x denotes "probability", taking into account that the process starts from x. Along these 
lines we have used the conformal invariance of the SLE measure, leading to invariance of probabil- 
ities under conformal transport. 

We now show that the right hand side of (|43p may be computed easily by establishing a relation 
to conditioned Brownian motion. Consider the process xt = \f^Bt of simple Brownian motion 
with diffusion constant k on the interval / = [— 7r, it], identifying its endpoints, which starts from 
xo = x £ I. We introduce the propagator via P(y,t; x)dy = P x [xt £ [y, y + dy]]. It is solution of 
the diffusion equation dP(y,t; x)/dt = (n/2)d 2 P(y, t; x)/dx 2 , and explicitly given by 

v n=— 00 x ' 

Let us condition xt to arrive at x p = ir (mod 2n) for some given time p > 0. Using elementary facts 
about conditional probabilities we see that the new process has a propagator Q(y, t; x) defined via 

Q(y, t; x)dy = P x [x t £[y,y + dy] \x p = ir (mod 2ir)\ 
P(ir,p-t;y) 



P(ir,p;x) 



P(y, t; x)dy, for < t < p. (45) 

The conditioning therefore leads to a drift [26( that can be read off from the diffusion equation for 
Q(y,t;x): 

dQ(y,t;x) d (( d \ \ Kd 2 Q(y,t;x) 

~ \ \ "a - haP(7T,p — t;y) Q(y, t;x) \ + — — (46) 



dt dy \\dy ' J J 2 dy 

Specialising to k = 2, we conclude that the conditioned process is solution of the stochastic differ- 
ential equation 

dx t = 2 (-^ In P(tt, p - t; x t )\ dt + V2 dB t 

= v(p-t,x t )dt + V2dB t (47) 

Hence we see that in the case k = 2 the motion Xt induced by the Loewner mapping gt(z) on 
the upper boundary of the cylinder is the same as Brownian motion xt on a circle, starting from 
x and conditioned to visit tt (mod 2ir) at time t = p. However, it is a simple exercise to compute 
the statistics of the winding number of the latter. Again, we shall say that xt has winding number 
N = n if its equivalent on the covering space M of / arrives at x + 2im at time t. From simple 
conditioning, we have 

exp (— (ir — x + 2Tin) 2 /2Kp) , , 

p, [Af = n|l, = 7T (mod 2,)] = EZ JJ p( _ {w _ - + ;;, i)V ' 2Krt (18) 

Summation over n from to 00 yields the announced result for k = 2 

u(px)-P\X >01- g = oexp(-(x- ff (2n + l))V4 P ) 



Thus we have shown that we may indeed reinterpret the winding of SLE2 in terms of the winding 
of a conditioned one-dimensional Brownian motion. In fact, this is the profound reason explaining 
why the preceding transformation leads to simplifications for k = 2. 



D. The winding of loop-erased random walks with fixed endpoints 

In this section we consider the case of a LERW with fixed endpoints from to x + ip, —tt < x < ir 
in the cylinder geometry. On the covering space, it is thus allowed to exit at x + 2ixN + ip with 
arbitrary winding number iV~ S Z. We shall be interested in the law of N which may be obtained 
by conditioning LERWs to exit at the given boundary points. 

Let us ask for the probability that the LERW has made N = n windings if we condition its trace 
to exit on a subinterval [a, b] + ip, —tt < a, b < tt on the upper boundary. It is given by 

p[*=»ir..» 6 m] = p[JV :" nRe ] p !,!'''''" = ■f> A ^+ 2 '"> , ,5„) 

1 P[Re>€ [a, 6]] /*diA 2 (p,x) 

where A 2 (p, x) and \ 2 (p,x) are the distributions defined in (|23l28p . Taking the limit a, b — > x 
amounts to forcing the LERW to exit at x + ip. It has been noted previously that (at least for 
simply-connected domains) this procedure of conditioning leads to chordal SLE2 from to x + ip 



271 . l28l . |29|. Supposing that this property holds for doubly-connected domains as well, we obtain 



the statistics of the winding number of SLE2 with fixed endpoints and x + ip on the boundary 
of T p : 

Po^WiV = n]= = ^ (51) 

2^ X 2 (p,x + 2nk) 

k=—oo 

where — tt < x < tt This result is non-trivial (at least for us), since there does not seem to 
exist an obvious way to recover this probability law from considerations of underlying random 
walks/Brownian motions. 

As an application of (|5ip we compute the generating function for the moments of N. For [i > 
let us write 



^ 00 
A 2 {p,x) ^ 



n=— 00 
1 d 



1 e w 2 /^ 2 q f e ^/ 27r e(p,x + pn/7v) ^ 
A 2 (p, x) 2 sinh(///2) dx \ e(p, x) 



(52) 



Strictly speaking, the calculation is valid for /j, > but the result extends to negative [/,, too. The 
result may be used in order to study the winding behaviour in the radial and dipolar limit. 



Radial limit. For p — > +00, we find a generating function for the cumulants given by 

Pfi 2 ( /Vg \ 

2tt 4vr 2 ° S Vsinh(^/2) J 



lo S <^""> = S + tl + log C-^r) + at,-'), (53) 



where we have used in that limit x) — ► l/(2-7r). As in section IlIIBI the cumulants (N n ) c with 
n > 2 converge to constants which can be obtained from the Taylor series expansion of the function 
log[/x/(2 sinh(/i/2))] with respect to (i. Corrections to this limit are exponentially small cx e~ p . 

Dipolar limit. Having in mind the discussion in section UlIBI we expect the probability distribution 
of N to be concentrated at N = as p — > + (independently of the value taken by |x| < ir), and 
the probabilities for events occurring with iV 7^ to be exponentially small oc e~ 2n ' p '. By means 
of explicit computation we find 

<e-^> = 1 + 4sinh (f) ( 2 sinh (^^) + sinh (^^)) (54) 
up to corrections of the order of 0(e _47r2//p ), what confirms the preceding discussion. 



IV. ENDPOINT PROBABILITIES FOR k = 4 ON THE CYLINDER 

It is surprising that we obtain the endpoint distribution on the cylinder k = 2 by periodisation of 
the equivalent dipolar result ([5]). For k = 4 we have applied the same idea to look for a solution 
in the form of a 2-7r-periodised result found in dipolar SLE4. These two values seem to be the only 
ones for which the periodisation yields exact results. The case k = 4 is related to the massless 
free field theory with central charge c = 1. In fact, SLE4 curves arise as zero lines of the field 
<p(z,z) emerging from discontinuities in Dirichlet boundary conditions. Moreover, SLE4 arises in 



the scaling limit of several lattice models such as the harmonic explorer and domino tilings 30|, l3l| . 
A variant of the harmonic explorer leading SLE4 in doubly-connected domains as discussed in this 
article was suggested in 0]. 

In this section, we shall show that the periodic endpoint distribution 



^ 1 



-°° cosh 



7r(x + 27T77.) 



2p 

is solution of (|12j) with n = 4. Let us notice, that (I55p corresponds to the endpoint distribution of 
two-dimensional Brownian motion on the cylinder with reflecting boundary conditions for y = 0, 
stopped for the first time when it reaches altitude y = p (see appendix IA 2\\ . 
We analytically continue A^p, z) to complex z£C and introduce the function 

<9A 4 (p, z) d d 2 A 4 (p,z) 
5(P ' Z) = dp' -^{v(p,z)A 4 (p,z))-2 — 2 (56) 

with v(z,p) as defined previously. If A^{p,z) is a solution of (|12|) with k = 4 then g(p,z) must 
vanish for z£l. We shall even show that g(p, z) vanishes for A^p, z) defined in (|55p for all z G C. 
A±{jp,z) is elliptic with periods u\ = 2tt and 0J2 = 4pi and antiperiodic with respect to the half- 
period W2/2 = 2pi. Thus, g(p,z) is elliptic for all p with the same periods u\ and U2 (as z — > 
z + Apt, the additional term arising from the quasi-periodicity of v(p,z) is compensated by the 
derivative of A with respect to p). We may therefore restrict our analysis to the rectangle defined 
via — 7r < Re z < tt and < Imz < 4p. The poles of A^PjZ), v(p,z) and their derivatives are 
located at z = ip and z = 3ip. The principal parts of the Laurent series expansions at z = ip for 
the different terms in (1561) read 



dA±(p,z) = 1 + d 2 A±(p,z) = 2i + 

dp tt{z — ip) 2 ' dz 2 ir(z — ip) 3 

d(v(p, z)A4(p, z)) 1 4i 



dz ir(z — ip) 2 tt(z — ip) 2 



+ 0(1) (57) 



Insertion into therefore shows that all divergent terms cancel out so that g(p, z) has a removable 
singularity at z = ip. The Laurent expansion around z = 3ip leads to the same result. In fact, we 
have 

dA A {p,z) = 3 + d 2 A,(p,z) = 2i + 

dp ir(z — 3ip) 2 ' dz 2 ir(z — 3ip) 3 

d(v (p, z)A4(p, z)) 3 4i 



dz tt(z — 3ip) 2 tt(z — 3ip) 3 



+ 0(1) (58) 



Therefore g(p, z) has a removable singularity at z = 3pi, too. Hence it must be a constant g(p, z) = 
A(p) what follows from Liouville's theorem However, because of g(p,z) = —g(p,z + 2pi) we 
find A(p) = for all p, so that g(p, z) = which is equivalent to say that for k = 4 the given form 
for A^(p, x) is the correct probability distribution function. 
Let us note that (|55|) may be written in terms of the Burgers potential e(p, x). 

A i (p,x)=M(p) e{P f~ 7r) . (59) 
e(j>, x) 

The normalisation factor M(p) may be found in a similar way as (|32H and (|33p . and is given by 

7re(p,ip-7r) dz it Enez e_n(n+1)p 

However, despite this very suggestive form of A^(p, x) as a ratio of two solutions to the simple dif- 
fusion equation we have not found algebraic simplifications in order to solve the partial differential 
equations for k = 4 with general boundary conditions. 



V. CONCLUSION 



In this paper, we have studied winding properties of loop-erased random walks around a finite 
cylinder by means of stochastic Loewner evolutions. Relating the computation of the endpoint 
distribution A(p, x) of SLE K on a cylinder to a problem of diffusion-advection of a passive scalar in 
a Burgers flow, we were able to explicitly determine \2(p, x) in the case k = 2. The behaviour in the 
limit of very thin and very large cylinders was studied and we pointed out a non-Gaussian behaviour 
of the winding properties for long cylinders. Furthermore, conditioning the loop-erased random 
walks to exit via a given boundary point, we were able to compute the probability distribution of 
the winding number for walks with fixed endpoints. We have shown that the somewhat surprising 
simplifications for k = 2 may be related to the winding of conditioned one-dimensional Brownian 
motion on a circle. Moreover, we have determined cylinder endpoint distribution in the case k = 4. 
A relation to reflected Brownian motion was pointed out. As for k = 2, it is related to the 
periodisation of the dipolar endpoint distribution. However, these seem to be the only values for 
k for which this property holds. It remains to see whether closed results can also be obtained for 
other values of k. We hope that our results can be useful to test recent conjectures 0] about SLE 
properties of interfaces in numerical studies in a cylinder geometry. 
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Appendix A: ENDPOINT DISTRIBUTIONS FOR TWO-DIMENSIONAL BROWNIAN 

MOTION 



In this appendix we compute the endpoint distribution for planar Brownian motion z% = xt + iyt 
on the cylinder, starting from z = and stopped as soon as it reaches the altitude y = p. For 
simplicity, we first consider the geometry of an infinite strip {z £ C|0 < lmz < p} and then 
periodise along the real axis in order to obtain the results on the cylinder T p . The diffusion in 
x— direction is unconstrained and has a propagator P Xo= o{x t £ [x, x + dx]] = exp(— x 2 /2t) dx/\/2-Kt. 
For the motion in y— direction, we consider absorbing and mixed boundary conditions. 



1. Absorbing boundary conditions 

Consider the motion in y-direction, starting from yo = e > 0. We shall impose absorbing boundary 
conditions at y = and y = p, and condition the process to exit at y = p. The propagator PA(y,t) 
defined via PA(y,t)dy = P yo=e [y t £ [y, y + dy]] is solution of the diffusion equation dPA(y,t)/dt = 
1/2 d 2 PA(y, t)/dy 2 and given by 

Pa(v> = - E sm {-j- ) 8111 {— ) ex P (—yr ) ( A1 ) 

We obtain the exit-time distribution /a(£) at y = p from the probability current at this point. 
However, if we condition the diffusion to exit at y = p, we furthermore must divide this current by 
the exit probability P yo=e [exit at y = p] = e/p. Hence 



fA(t) 



P dP(y,t;e) vr ^ +1 . f nire\ / ir 2 n 2 t 

= — > -1 nsm exp t— 

2e dy y=p pe ^ } \ p J P \ 2p 2 

9 00 / —2,^2 



2 (JXJ / 2 2j.\ 

^D-D^exp , as^0+. (A2) 

n=l 

We obtain the distribution ^a{Pi x ) of the endpoint by integration of the free propagator in in- 
direction weighted by /a(*) with respect to time t: 

f°° dt ( x 2 \ 

OO _ 1 

p^ K ' 4p(cosh (irx2p) 2 v ; 

k=l 

Periodisation in x with period 2tt yields the endpoint distribution on the cylinder A2(p,x) (I23p . 



2. Mixed boundary conditions 



Consider the same problem as above, but impose reflecting boundary conditions at y = and 
absorbing boundary conditions at y = p. Let PM(y,t) the propagator in y-direction, defined 
via PM(y,t)dy = P yo=0 [y t £ [y,y + dy]]. It is solution of the diffusion equation dPM{y,t) / dt = 
1/2 d 2 PM(y,t) /dy 2 and given by 

„ / \ 2^ /(2n+l)vry\ / (2n + l) 2 vr 2 A , A 

PM(t,y) = -gooB (^-^J exp (- 1 J j • (A4) 



As previously, we compute the exit-time distribution from the probability current at y = p: 

00 / CO„ 1 i\2^2 



fMit) = £ J2(- i r(2n + 1) exp (- (2n+ J 2 )7rt ) • (A5) 



n=0 

Consequently the distribution of the exit point at the upper boundary of the strip reads 

f°° dt ( x 2 \ 

mm(p ' x) = / 7^rt exp {- Yt) fM{t) 

1^, / (2n+l)7r|x|\ 11 , K . 

= ->(-! exp [ — i — ; 1 1 = ; — A6) 

P^ V 2 P J 2pcosh(irx/2p) y ' 

Finally, periodisation with respect to x with period 2tt leads to A^{p, x) for the cylinder T p (|55l) . 



Appendix B: ENDPOINT DISTRIBUTION FOR HEXAGONAL LATTICES 



The purpose of this appendix is to give some idea about finite-size corrections to the scaling limit on 
the cylinder. Here we evaluate the discrete equivalent to f|23|) and compute the leading correction 
due to lattice effects. In principle, we dispose of a great variety for the choice of the lattice. Finite- 
size corrections not only turn out to be dependent on the choice of geometry, but on the lattice 
type as well. We shall content ourselves with the honeycomb lattice - the most common choice. 




Figure 4: (a) Cylindrical lattice domain of 5 x 3 honeycombs with walk from lower to upper boundary. The 
highlighted curve represent some walk from the lower to the upper boundary (b) Coordinate system for 
equivalent "brick wall" lattice. 



We consider a finite lattice tube of m x I honeycombs, as shown in figured^. We decompose the 
honeycomb lattice into two triangular sub-lattices whose sites are coloured by • and o. Nearest 
neighbours of a •— site thus belong to the o— lattice and vice versa. We study lattice walks starting 
from a •— site at the lower boundary to some any o— site at the upper boundary and shall impose 
absorbing boundary conditions on o— sites at either boundary. In fact, it shall be convenient to map 
the problem to a "brick wall" lattice as shown on figure Hb, and introduce a suitable coordinate 
system. Then it is straightforward to write the the master equations for the probability P./ (x, y; t) 
that the walker starting from (0, 0) can be found at (x, y) at time t: 

P (x,y;t + 1) = 1 (P.(x + l,y;t)+P.(x - l,y;t) + P.(x,y - l;t)) (Bl) 

P. (x, y; t + 1) = - (P (x + 1, y; t) + P a (x - 1, y; t) + P Q (x, y + 1; t)) (B2) 



with boundary conditions P Q (x,y = 0;t) = P (x,y = l;t) = and P,(x,y;t = 0) = 5 X ',o^,o where 
8ij denotes the Kronecker symbol which takes the value 1 if i = j, and otherwise. Moreover, 
identification of x and x + 2m leads to periodicity P Q (x + 2m, y; t) = P Q (x, y; t), P,{x + 2m, y; t) = 
P,(x,y;t). We shall need the local occupation times G,/ Q (x,y) = Ylt^o P»/o(%, V'-, £)■ The lattice 
walker is stopped if it walks from some •— site at y e = 1 or y e = I — 1 to a o— site at y e = 
or y e = I respectively, what happens with probability 1/3. Therefore, the exit probability at the 
upper boundary is given by 

P[x e = x,y e = l] = ~G.(x,l-l), (B3) 

where x e denotes the (random) horizontal position of the exiting walker. Notice that x e only takes 
odd/even integer values for even/odd /. Combining (|B1| IB2[) . it is possible to eliminate G (x,y) 
and find the the second-order difference equation 

G.(x,y) =^(G.(x + 2, y) + G.(x + l,y + 1) + G.(x + l,y - 1) 

G.{x - 2,y) + G.(x - l,y + 1) + G.(x - l,y - 1)), (B4) 

for < y < I — 1. Using the boundary conditions, we furthermore find the equations 

G.(x, 0) = l (G.(x + 1,1)+ G.{x - 1, 1)) +1 4,o (B5) 

o 

G.(x,/-1) = i(G.(x + l,/-2) + G.(x + 2,/-l) 

+ G.(x-l,l-2) + G.(x-2,l-l)) (B6) 

Solving (|B4p . (|B5j) and (|B6jl all together leads to a rather lengthy and tedious calculation. The 
reader may find a detailed account of the strategy with applications to other lattice types in [33]. 
Here we only sketch the solution. 

First, we solve the bulk equation (|B4[) by a separation ansatz G,(x, y) = P(x)Q(y). Forming linear 
combinations of solutions found from this ansatz, we obtain the general bulk solution 

G.(x,y) =A + C y + (A m + C m y)(-l) x+ y (B7) 

+ (A k e iakX + B k e- iakX )(e^ y + C k e~™) 
keh 

+ J2 {-l) y (A k e iakX + B k e~ iakX ){e lky + C k e~^ ky ) 

where Ai, Bi and C{ are constants, I\ = {k < k < m/2 V 3m/2 < k < 2m} and I2 = {k G 

Z I m/2 < k < 3m/2, k 7^ m/2} index sets, and 

irk 3 — cos 2a k 

a k = — , cosh7fc = — j j— and k = 0, 1, . . . , 2m, — 1 (B8) 

m 2|cosa/c| 

However, we exclude k = m. Moreover, for some technical reason the formula only is valid as long 
as m is odd. The remaining constants are determined by matching (jB7[) to ()B5p and (jB6[) what 
yields 

1 ^_ _ „27fc(l-i) l + 2e^|cos« fc | 



^4fc — -Bfc — 7; — i 1 . 1 / — 7T, Cfc — — e 



2m| cos sinh(7fc/) ' 1 + 2e | cos | ' 



and Aq = A m = (31 — l)/(4mZ), Co = C m = —3/ (4ml). Putting all pieces together, we find the 
solution for (|B3p . However, if we condition the walk to exit at the upper boundary, we still have 
to divide by the probability P[y e = I] what amounts to a normalisation. After some algebra we 
obtain the final result (for odd m) 

(m-l)/2 

P[x e = x\y e = J]=i+ Y 21 c ^i^pK . (B9) 
m msinh7i.t 

Having this lattice result, it is interesting to compute its scaling limit and study corrections. In 
the case where m is odd, x takes even values from — (m — 1) to m— 1. Introducing a lattice scale 
a such that x = £/a, £ G R, we would like to study the scaling limit a — > + , m,l —> +oo in such 
a way that the geometrical height p = y3(l — l/3)a and width 2tt = 2am remain finite. Therefore 
l/m = p/ir — a I (\p3ix). 

1 y/3l ^> k cos/c£ 

X(p,t,a) = — P\x e = x\y e = I \ = > 1= 

V ' 2a L ltfe J nm ^^8^(^/3(1 /m)k) 

1 ^ i 



4V3(i/m) ^ Uhfi^I 
L V 2v^3(i/m) 



l + ^=-J A 2 (p,0+O(a 2 ) (BIO) 

We see that the first correction to the scaling limit is completely naturally proportional to the 
lattice scale a. 



Appendix C: ANALYSIS IN THE BULK 

This appendix collects some facts about bulk properties for SLE on the cylinder. Let us denote 
the probability, that the trace 7 passes to the left of some given point z on the covering space by 
u>(p,z,z). It is solution of the diffusion equation 

du(p,z,z) duj(p,z,z) du(p,z,z) n ( d d \ 2 

-H(p,z) + H(p,z) — + ^(piZ^) (CI) 



dp dz dz 2 \dz dz 

with boundary conditions u>(p,z,z) — » as Rez — > —00 and u(p,z,z) — ► 1 as Kez — » +00. 
The derivation of ()Cip is similar to (jlip and (|12p . using an infinitesimal argument. Using the 
relationship H(p, z) = v(p, z — ip) — i to the velocity field defined above and (fT6j) . it is not difficult 
to show that H (p, z) is a Burgers flow, too. We have 

9H(p,z) dH(p,z) d*H(p,z) 

dp ~ H ^ Z > Qz + dz* (L2) 

Using the Cole-Hopf transformation (|18|) . we may write H(p,z) = dU(p, z)/dz with U(p,z) = 
2 log e(p, z — ip) — iz — p/2. 

For k = 2 as Im z —* p~ we require that z, z) tends to A 2 (p, x) (with x = Re z), and u;(p, z, z) — > 
0(x) as Imz — » + . Although we have not found the solution with these boundary conditions, we 
have determined a special solution iv(p, z,z) with a periodised version of these boundary conditions 
which has interesting relations to dipolar SLE2. We ask that for Imz — > p~ it becomes &(p, z,z) = 
Jq dxA2(p,x) = (pv(p,x) +x)/2tt, and that for Imz — » + it reproduces the quasi-periodic step 
function uj(p, z,z) = lim^oo J]fe=-n( si g n ( x ~ 27rfc))/2. 



Obviously, both H(p,z) and H(p,~z) are solutions to (jCip . however with a simple pole at z = 0. 
Inspired from the boundary condition at Imz = p, we shall consider the function f(p,z,~z) = 
(hnz)(H(p, z) + H(p,z)), which has no pole. However, using the differential operator A defined as 



A 



dp ^' Z ^ dz 



H(p, z] 



d_ 

dz 



8_ d_ 

dz dz 



(C3) 



we see that Af(p, z,z) = — (H(z) 2 — H(z) 2 )/2i. Hence we should find some solution to Ag(p, z) = 
iH(z) 2 /2 without poles (since this would again lead to a singular solution). Hence we have to 
search for a suitable g(p, z) with at most a logarithmic singularity at z = 0. This can be done by 
observing that U(p, z) = 2 log e(p, z — ip) — iz — p/2 is solution to the equation 



dll(p, z) 
dp 



1 f dujp, 

2 V dz 



+ 



d 2 U(p, z) 
dz 2 



1 



, dU(p,z) d 2 U{p,z) , 
H{p > z) ^z- + ^ 2 —-2 H{P > z) 



(C4) 



as can be seen from Burgers equation for H(p, z) (lC2j) . Hence we have found a solution uj(p, z,z) = 
C\ + C2(f(p,z,z) — i(U(p,z) — U(p,z))) 1 where C±, C2 are constants. It is a solution of (jCID 
without poles for < Imz < p. The choice C\ = 1/2 and C2 = 1/4tt leads to the desired boundary 
conditions for u)(p,z,z): 



I Im z % 

u>(p,z,z) = - + —(H(p,z)+H{p,z)) + —(U{p,z) -U(p,z)) 



2 4vr 

u(p, z,z) is related to ui(p,z,z) by periodisation 



4tt 



n , 1 

u(p,z,z)= lim N \ uj(p, z + 2nk, z + 2irk) 



(C5) 



(C6) 



i.e. we have only found a (quasi-)periodic version of uj(p, z, z). 

The dipolar limit p — > + turns out to be quite interesting. In that limit the quasi-periodicity is 
irrelevant and one obtains the result for dipolar SLE2. One can replace in that limit H(p,z) — > 
H ss (p,z) = ^cothirz/2p — z/p, the single shock expression, and the corresponding expression 
21oge(p, z) ~ —z 2 j{2p) + 21ogcosh7rz/2p, up to a p-dependent unimportant constant. One then 
easily obtains from f|05|) : 



1 TTZ Imz ( TTZ TTZ 

1 Im log sinh 1 coth h coth — 

vr B 2p Ap \ 2p 2p 



(C7) 



This is nothing but the probability that dipolar SLE2 passes to the left of a given point z within 
a strip of height p. It satisfies the general equation: 



7T ^ TTZ d IT TTZ d K ( d d 

2p dz p 2p dz 2 \dz dz 



P 



"dip, z, z) = 



(C8) 



studied in [13J where solutions could be found only for k > 4 (in the form of simple harmonic 
functions). Here we recover, from a limit of a more general object defined on T p , a result for 
k = 2 obtained only very recently by a rather different method 34|. The fact that (|C7p satisfies 
both (|C8j) (an equation with no d p terms) and the dipolar (i.e. single shock) limit of (|Clj) (an 
evolution equation as a function of p) is easily understood by noting that it also satisfies (pd p + 



zd z + zdz)oJd{Pi z i z) = which expresses the dilatation invariance of the dipolar SLE, explicitly 
broken by the period of the cylinder. 

One may wonder whether there is a way to extend the known harmonic solutions [l3| of the 
dipolar limit, to the full cylinder. Here we only present the case k = 4 which can be constructed 
by integration of the result given in section HVl In fact, this amounts to periodise the result from 
dipolar SLE for a strip of height p 

1 7TZ 

Udip, z, z) = 1 Im log tanh— (C9) 

in an appropriate way. Symmetric periodisation and leads to 

Q(p,z,z)= lim Y (I --Im log tanh Mi+^E*)^ (C10 ) 

which is a harmonic solution to (|Clj) with required boundary conditions at z = x + ip 

u(p,x + ip,x — ip) = / dxA4(p,x). (Cll) 
Jo 



Appendix D: PATH INTEGRAL FOR GENERAL k 

Let us indicate how one can take advantage of ([24"]) to express the solution of (fT2j) as a path integral 
for general k. We may write 

/oo 
dy K(p,x;y)ip(0,y) (Dl) 
-oo 

where K(p,x;y) is the Euclidean propagator for (|24p . It given by a path-integral representation 

f'X {jp ^ — X 

K(p,x;y)= ' [dx(t)]e~ s[xit)] (D2) 

Jx{0)=y 

where the action is defined as 

sm]= !\ t m + ^^m. (D3) 







2k 2k dx 

Hence, we find a nice path-integral representation [la . for the general solution uj(p, x) of ([15 



dye(0,y) 2 /^(0,y) / ^ [6x{t)]e~^ (D4) 

-oo J x(0)=y 

which reproduces the exact solution (|27p in the case k = 2. It may be useful for perturbation 
theory around k = 2, using the Fourier series decomposition: 

«^ 4 e ~ np 

v(p, x) = 2_^ - _ e _ 2np sin(nx) (D5) 
n=l 

if one treats the apparent singular behaviour of e(0, y) 2 ^ K for k / 2 (e.g. shifting the integration 
from to small po an d using the known dipolar result for u>(po,x)). 
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